1 



Quenched charmonium spectrum on anisotropic lattices * 

CP-PACS Collaboration : A. Ali Khan? S. Aoki b R. Burkhalter a - b S. Ejirif M. Fukugita^ 

S. Hashimoto? N. Ishizuka a ' b Y. Iwasakif- b K. Kanaya b T. Kanekof Y. Kuramashi? K.-I. Nagaif 

M. Okamoto b M. Okawaf H.P. Shanahanft Y. Taniguchi b A. Ukawa a ' b and T. Yoshie a - b 

a Center for Computational Physics, University of Tsukuba, Tsukuba, Ibaraki 305-8577, Japan 

b Institute of Physics, University of Tsukuba, Tsukuba, Ibaraki 305-8571, Japan 

Tnstitute for Cosmic Ray Research, University of Tokyo, Kashiwa 277-8582, Japan 

d High Energy Accelerator Research Organization (KEK), Tsukuba, Ibaraki 305-0801, Japan 

We present the results of quenched charmonium spectrum for S- and P-states, obtained by a relativistic heavy 
quark method on anisotropic lattices. Simulations are carried out using the standard plaquette gauge action 
and a meanfield- improved clover quark action at a^ 1 = 3-6 GeV with the renormalized anisotropy fixed to 
£ = a s /a t = 3. We study the scaling of our fine and hyperfine mass splittings, and compare with previous results. 



1. Introduction 

Conventional formalism of lattice QCD fails to 
describe heavy quarks because auiQ >0(1) for 
charm and bottom quarks on current lattices with 
a -1 ~ 1-3 GeV. In heavy hadrons, however, typ- 
ical scale for the spatial momenta are compara- 
ble with that of light hadrons. A momentum 
with O(toq) appears only in the time direction. 
Therefore, as proposed by Klassen, it is natural 
to adopt an anisotropic lattice with a small tem- 
poral lattice spacing with atrriQ <C 1 for a lattice 
QCD description of heavy hadrons jj]|2| . 

In this report, we present results of our on- 
going calculation on the quenched charmonium 
spectrum using the anisotropic method, which at- 
tempts to test the feasibility of this approach. 



Table 1 

Simulation parameters. Lattice spacing is fixed 
by the Sommer scale ro = 0.5 fm. 
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2. Method and simulation parameters 

The gauge action we use is given by 

where £o is the bare anisotropy. The renormal- 
ized anisotropy, defined by £ = a s /at, can be de- 
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termined by Wilson loops. In this study, we use 
the result of § for £(&,/?) • 

For the quark action, we employ an O(a)- 
improvcd Wilson-type quark action on an 
anisotropic lattice QJ^Q: 

Sf = ^2{4>xipx 
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Table 2 

S- and P-state operators. 

state J^ u name local non-local 

3 S 1 1— J/i> fy^ 
1 Pi l+~ h c ipajjip i/ry 5 Aiip 
3 Po 0++ Xco H> tY.ilM 
3 Pi 1++ Xd ipli-ysip fHAj - 7jAj}V 
3 P 2 2++ X c2 HliAi - (E) 
^{7A-+7jAi}^ (T) 



Figure 1. The speed of light and effective speed 
of light for r\ c at tuned £ and /3 = 5.9. 
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where i4T s t and c S; t are the spatial and temporal 
hopping parameters and clover coefficients. The 
bare quark mass is given by 

a t m q0 = l/2K t - 3/C - 1, (3) 

where ( = K t /K s . We adopt the meanfield- 
improved clover coefficients: 

f l + £ 1 
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where (U s ) — (Pss') 1 ^ 4 with P ss i the spatial pla- 
quette and (U t ) = 1. 

Our simulation parameters are summarized in 
Table |l|. We study lattices with a fixed renor- 
malized anisotropy £ = 3 and approximately the 
same spatial lattice size of about 1.6 fm. In order 
to study the charm quark, we simulate three lat- 
tices with the lattice cutoff in the range a^ 1 = 3-6 
GeV. 

The parameter £ has to be tuned such that the 
relativistic errors are eliminated. For this purpose 



preparatory simulations are performed at ( = 2.8, 
3.0 and 3.2. We measure the pseudoscalar meson 
mass at four lowest on- and off-axis momenta and 
perform a fit of form, 



E(p) 2 = E(0) 2 + c 2 (p)p 2 



(5) 



to extract the "speed of light" , c(p) , using three 
or four lowest momenta. The tuning condition for 
£ is c(p = 0) = 1 which requires the relativistic 
dispersion relation to be restored for small mo- 
menta. (For massless free fermions, c(p = 0) = 1 
is satisfied at C — £•) Results for ( are summa- 
rized in Table [l| 

Figure [l] shows a typical result for the effective 
speed of light c cS (p) = ^E(p) 2 - E(0) 2 /p, ob- 
tained by a simulation made with a tuned value 
of £. The wide plateau indicates that the linear- 
ity of E{p) 2 in p 2 is well satisfied. The horizontal 
solid line is the result of c(0) from a global fit 
according to (j|), where the dashed lines indicate 
its error. We see that c(0) is consistent with 1. 
As summarized in Table [j], the condition c(0) = 1 
is confirmed to be satisfied within 1-2% with our 
values of £. 

2.1. Measurements 

At each /?, hadronic measurements are made on 
configurations separated by 100-200 iterations, at 
two values of m q o . See Table [l]. These two quark 
masses are chosen such that the charm quark 
point can be interpolated. 

We measure all S- and P-state mesons using 
both local and non-local operators compiled in 
Table [|. For the non-local operators we adopt an 
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Figure 2. The charmonium mass spectrum at £ = 
3 and (3 = 5.9. The scale is fixed through the 
Sommer scale rrj. 

Table 3 



Preliminary continuum estimates of mass split- 
tings and the experimental values. 
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exponentially smeared derivative sources: 

f % (x) = A s e- B ^-A s e- B ^ (i = 1,2,3),(6) 

where A s and B s are the smearing parameters. 
We extract meson masses by a single cosh fit. Er- 
rors are determined by a jack-knife method. 

3. Charmonium spectrum 

We fix the scale using either the Sommer scale 
ro = 0.5 fm or the l 1 Pi — IS splitting. The charm 
quark point is then defined by setting the spin 
averaged IS meson mass (IS). Figure || is a typi- 
cal result for the charmonium mass spectrum ob- 
tained at (3 = 5.9. In the following, we study 
the lattice spacing dependence of the spectrum. 
With the 0(a) improved action, the leading scal- 
ing violation is expected to be proportional to a 2 . 
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Figure 3. The scaling behavior of hyperfine- 
splitting AMjw.^. Filled circle and square are 
our results with scale from the Sommer scale ro 
and the IP — IS splitting respectively. Open di- 
amond and triangle are results of Ref . 0] and || . 
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Figure 4. The scaling behavior of AM Xcl _ Xc0 . 

3.1. Hyperfine splitting 

Figure |3] is our result for the S-state hyperfine- 
splitting J/ip — t] c as a function of the lattice 
spacing. Results using the scales from ro and 
l 1 Pi — IS are presented. We find that, although 
our results for the hyperfine-splitting can be ap- 
proximately fitted by a linear function of a^, the 
two results from different scales lead to different 
values even in the continuum limit. This discrep- 
ancy should be regarded as a quenching artifact. 

Plotted together are the previous lattice results 
by Klassen || and Chen|5) at £ = 2 and 3, using 
the same lattice action. A difference between our 
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3.2. Fine structure 

Results for the P-state fine structure Xci — XcQ 
are shown in Fig. |J. Although the scaling behav- 
ior of our results is similar to those of and ||, 
the continuum limit deviates by about 2a for the 
data using the ro scale. 

We also show the results for the Xc2 — Xci split- 
ting and the ratio \c2 - Xci/ Xci - Xco in Figs. | 
and g. No data from other groups are available 
for these quantities. 
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Figure 6. The ratio AM Xc2 _ Xcl /AM Xcl 
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simulation and those in || and || is the choice of 
the tadpole factor for the clover coefficients; we 
use the fourth root of the plaquette expectation 
value, while the mean link in the Landau gauge is 
used in Q and ||. To check the consistency, we 
also performed a simulation at /3 = 5.7 using the 
Landau gauge mean link for the meanfield. As 
shown in Fig. |^ by a filled diamond, our result is 
consistent with those by Klassen. 

We find that our result from the scale ro is 
approximately parallel to the other two results, 
and lead to a continuum limit which is about 10% 
lower than that of the others. The origin of this 
discrepancy is not clear to us at present. 



4. Conclusions 

We have presented our results for the charmo- 
nium mass spectrum obtained on an anisotropic 
lattice with £ = 3. Our results for fine and hy- 
perfine mass splittings are shown in Figs.|||| Our 
preliminary values from a continuum extrapola- 
tion linear in a 2 are summarized in Table 3. We 
find that these results arc much smaller than ex- 
perimental values. 

We also find that the results are quite sensitive 
to the choice of the clover coefficient. A naive 
continuum extrapolation, linear in a 2 , does not 
resolve the discrepancy between different choices 
of the clover coefficient. We are now extending 
the simulation to a finer lattice, to examine the 
reliability of continuum extrapolations. 
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